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The one loop effective action in a Schwarzschild bacfcground is here used to compute the Zero 
Point Energy (ZPE) which is compared to the same one generated by a gravastar. We find that only 
when we set up a difference between ZPE in these different background we can have an indication on 
which configuration is favored. Such a ZPE difference represents the Casimir energy. It is shown that 
the expression of the ZPE is equivalent to the one computed by means of a variational approach. To 
handle with ZPE divergences, we use the zeta function regularization. A renormalization procedure 
to remove the infinities together with a renormalization group equation is introduced. We find that 
the final configuration is dependent on the ratio between the radius of the wormhole augmented by 
the "brick wall" and the radius of the gravastar. 
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; I. INTRODUCTION 

, Black holes are well accepted astrophysical objects by scientific community. The simplest example of a black hole 
is the spherically symmetric vacuum solution of the Einstein field equations, 

o^: ds' = " " ^) '^^^ + 1 - ^Lg + + ^'^'^') ' (1) 

known as the Schwarzschild solution. Despite of many theoretical and observational successes, a number of paradoxical 
problems connected to black holes also exist|l)], which frequently motivate authors to look for other alternatives, in 
which the endpoints of gravitational collapse are massive stars without horizons. In 2001, Mazur and MottolaQ] 
proposed an alternative model to black hole as a different final state of a gravitational collapse. In this model, the 
strong gravitational forces induce a vacuum rearrangement in such a way to avoid a classical event horizon. A phase 
transition is associated to the quantum gravitational vacuum together with a topology change. Such a model has 
. been termed gravastar (gravitational wacuum star) and it consists of three different regions with three different Eqs. 
of state 
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, /. Interior : < r < ri 



II. Thin Shell : ri < r < ra p = +P, (2) 
///. Exterior : 7'2 < r p = p = 0. 

At the interfaces r = ri and r = r2, we require the metric coefficients to be continuous, although the associated first 
derivatives must be discontinuous. Globally the metric can be cast into a form very close to the Schwarzschild line 



' element ([T]) 



where 



dr"^ 

ds'^ = (r) dt^ + ^ + {d0^ + sin^ edc/)^) , (3) 



N^r) = l-^r^ and b{r) = ^r^ (4) 



for the interior region of a static de Sitter metric, while 
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(r) = 1 and b (r) = 2MG (5) 



'Electronic address; |Remo.Garattini@unibgTTt] 



2 



for the exterior region described by a Schwarzschild spacetime. The intermediate region is represented by a thin shell 
endowed with a Minkowski metric. However, the thin shell is by no means necessary to obtain a gravastar. Indeed, 
DeBenedictis et al. Q proved that the shell region can be eliminated and the de Sitter spacetime can be directly 
joined to the Schwarzschild metric. This picture was also considered by Dymnikova 4] without invoking the term 
gravastar" . From the asymptotic point of view, a gravastar and a black hole share the same Arnowitt-Deser-Misner 
(AVAi) massQ, therefore at very large distances they appear the same object to an observer. One may wonder how 
can we distinguish a black hole from a gravastar. Different proposal have been considered. Harko et al.[6] suggested 
to compare the thermodynamical and electromagnetic properties of the accretion disk around slowly rotating black 
holes and gravastars. Another proposal comes from Chirenti and RezzollaQ, where the use of axial perturbations 
and the analysis of Quasi Normal Modes seems to show how to distinguish a black hole from a gravastar. Cardoso et 
al.Q discussed the ergoregion instability produced by rapidly spinning compact objects such as gravastars and boson 
stars with the results that ultra-compact objects with large rotation are black holes. On the other side, Chirenti 
and RezzollaQ found that stable models can be constructed also with J/M^ ~ 1, where J and M are the angular 
momentum and mass of the gravastar, respectively. Note that the comparison between a gravastar and a black hole is 
at the classical level, without any quantum contribution which should also be the source of the desired phase transition 
to make a gravastar. Note also that, in principle there exists another source of ambiguity: indeed besides a gravastar 
and a black hole, the exterior metric can be associated to a wormhole. In this paper, we would like to compare a 
gravastar to a wormhole. From the asymptotic point of view even the wormhole shares the same ADM mass. Hence, 
it should be very important the comparison of such objects invoking the Zero Point Energy (ZPE) contribution. 
However, it is well known that every form of ZPE also contributes to an induced cosmological constant. If we indicate 
with Eq^ the gravastar ZPE and with , the wormhole ZPE, in principle one can discuss the following inequalities 

«-i?o«^)|0 (6) 

establishing which geometry is energetically favored compared to the other one. Essentially, this is a Casimir-like 
calculation and lower the ZPE, the more stable is the final configuration. This inequality is also related to the decay 
probability per unit volume and time F, which is defined as 

F = Aexp i-Ici) = exp {-Ig [g^,]) J Vh^, exp (-if^ [/i,,,]) , (7) 

for an Euclidean time. Indeed, at least a tree level, we find that 

F ~ exp {-1,1 ) = exp (- {E§^ - E^) At) , (8) 
where Ar is the Euclidean time interval. Basically, the first step is the evaluation of the following expectation value 

provided one can compute the effective action related to Sg. To do calculations in practice, we fix our attention to 
the standard Einstein action without matter fields and with a cosmological term 



S = Sg + Sec. (10) 



where 



The least action principle leads to the Einstein's fields equations with a cosmological term. Nevertheless, never forbids 
to consider the cosmological term as the desired induced quantity by ZPE. Therefore, Eq.(|9l) must be modified into 

If we define the path integral 

^ = y 2? [g^iv] exp iSg [g,,^] (13) 
and we consider a gravitational field of the form 

g^iu = gp.u + hi_,^, (14) 
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then Eq. p^ becomes 

J Vg^i, expiSg [g^^] = expiSg [g^,.] Z2 = cxpiSg [g^^,] J 'Dhf_,^ expiSf^ [h^,^] , 
where we have assumed that the background (jfj^i, is a solution of the Einstein field equations and 
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with K = SttG. O^f""^ is a symmetric tensor operator with 



Sht^PSh'"" 



(15) 



(16) 



(17) 



and is the quantum fluctuation with respect to the background gau- After some integration by parts, Eq.([T6 
becomes 



5(2) = J_ 
2k 



M 
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1 



-/i^- (Az/i)^, + -h^h + h^PRp.h; - -h^''h%^,. 



R 



(18) 



Rfiah^ -\- Ri/ah^ 



A ^stands for the Lichnerowicz operator defined by 

(Al/i),,, = -V"Va/ip. - 2Rp^,0h''^ + R^o^K + ^.a/i;: = A/i^, - 2i?^„./3/i"'' 
where we have introduced the positive definite differential operator 

A - -V°Va. 

S'^^-' simplifies considerably when we are on shell, namely R^p — 0. Nevertheless, for future purposes, it is convenient 
keeping such terms in the expression of the Lichnerowicz operator and in Ea. (|18p . To extract physical informations 



(19) 



(20) 



from expression (|15p . we need an orthogonal decomposition which is equivalent to the Faddeev-Popov procedure, at 
least to one loop. From Appendix!^ we obtain[l3, 



Ti-ioop = ^ [rr In Af- Tr In A^/i] = ^ y" d'^xy/^ J 



d'^k 



(27r; 



4 In ^TT ^ 



d^k 



{2n\ 



4 In Ay J 



(21) 



where X'^rp and \^±_ are the eigenvalues of the Lichnerowicz operator for TT tensors and the transverse vector operator 
respectively. With the help of Eq. ((2T|) , Eq. (fT2)) simply becomes 



and if we identify 



A 

8^ 



4 m Arpj' — 



(27r 



(2^ 



4 In Ayj 



d'^k 
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TT 



d'^k 
(^ 



InAf 



(22) 



(23) 



we can interpret the cosmological constant as induced by quantum fluctuations of the gravitation field itself. Therefore, 
it is clear that inequality © can be directly measured by the induced cosmological quantity of Eq. . Moreover this 
identification will be useful for the removal of divergences. A first observation about Ea. (|21l) is in order. Note that 
nothing has been said regarding the famous conformal factor problem. From this point of view, we adopt the approach 
of Mazur and Mottola 12| in decomposing the gravitational perturbation. The super-metric free parameter "C" of 
Eg. dASI) leaves us the freedom to select the correct range in such a way the functional integration be convergent. 
Coming back to the Lichnerowicz operator (jl9p . one immediately recognize that finding the eigenvalues is not a trivial 
task in general. It is therefore necessary to adopt a convenient choice to manage Eg. (1221) . In a previous work we 
approached the cosmological constant problem with the help of the Wheeler-DeWitt Equation cast in the form of 
a Sturm-Liouville problem [isj. Essentially the cosmological constant is reinterpreted as an eigenvalue, calculated in 
a Hamiltonian formalism breaking the covariance of space-time. In the next section, we adopt the same strategy 
provided one looks at the true degrees of freedom. The paper is organized as follows: in section [TTl we reduce the 
effective action by restricting the modes of the perturbation, in section lTIIl we evaluate the functional determinants by 
means of a W.K.B. method, in section lTVl we compute the ZPE energy for the gravastar and the wormhole respectively 
and we compare them. Finally, in section |V] we conclude. 
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II. REDUCING THE ONE LOOP EFFECTIVE ACTION IN 3+1 DIMENSIONS 



How the Lichnerowicz operator decomposes in 3+1 dimensions it depends on the way one separates space from 
time. The AT>M variables offer a vahd example of such a decomposition. In terms of these variables, the metric 
background written in Eq.(I3]) becomes 



-N^dt^ + g,j {N'dt + dx') {N^dt + dx^) 



(24) 



We recognize that the lapse function N is invariant and the shift function Ni is absent. To have an effective reduction 
of the modes, we consider perturbations of the gravitational field on the hypersurface Y, <Z Ai. This means that we 
are freezing" the perturbation of the lapse and the shift functions respectively. In summary. 



N - 
N,, 



N 
> 



(25) 



corresponding to a restriction of the modes we are looking at. Choice 

hof, ^ (^ = 0, . . . , 3) 

in Eq. p^ . which can be reduced to 



is equivalent to set 



(26) 



5(2) = 1- 
2k 



d*Xy/—g 



M 



(27) 



Note that the modes we have eliminated satisfy the transverse traceless condition. The remaining modes are described 
only by spatial indices which are raised and lowered using gij and g*^ . Christoffel symbols and Riemann tensor are 
entirely constructed with the help of the three dimensional background metric. It is clear that even decomposition 
(jAip is affected by the reduction ([25)) which induces a rearrangement of the Ea. (P7)l . After a lengthy algebraic 
manipulation we arrive at 



^(2 



dt / d-^xNy/g 



-a An 



(28) 



where 



A^^ 



{Adi' 



(29) 



and 



Act = At 



12 2iV2 Q^2 



a. 



(30) 



It is immediate to recognize that for Einstein background Rij — Agij, cross terms vanish. Unfortunately, the 
Schwarzschild metric in three dimensions does not fall in this case. Nevertheless, the linearized action can be repre- 
sented in a short way on a suitable tensor space 



Q[{ij,m){kLn)] 



(31) 



with o[(»i.'")(fc'.")] a [(3 x 3) + 1] x [(3 x 3) + l]-matrix differential operator whose first (3 x 3) block matrix act on 
transverse traceless spin two field /i^ and whose last columns acts on the spin zero field a. The corresponding matrix 
can be read off from 







3 . - 



(32) 



To write the corresponding functional determinant, we observe that the following relations are valid for arbitrary 
triangular matrix operator[l4|: 



In det 



A C 
B 



= Tr In 



A C 
B 



= Tr 



InA C 
InS 
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^\nA + \nB = In dot 



Thus, the mixing term does not come into play and we get 



A 
B 



(det,.,, [o(*'^")]) 



= det I A, 



det^ A 



The same problem appears for the Jacobian. In Appendix I A 11 we show that Eq. (PT|) reduces to 



if f d k 



M 



and Eg. ip^ changes into 



■InA^ 



{T^i:' ')=~7^9i.v J ^^lnA^y + ^(5^^ + u^u^) j -^^...^y^, 



(33) 



(34) 



(35) 



(36) 



where is a time-like unit vector. This means that the energy density in 3+1 dimensions is not affected by the 
vector part which contributes only on the pressure terms^. Eq.([36]). 



III. W.K.B. APPROXIMATION OF THE FUNCTIONAL DETERMINANTS 



To evaluate pes and pw, we extract the energy density from Eq. (p6| and we get 



A 
8^ 



d'^k 



(2^) 



4 In ^TT 



d^k f+°°duj^ / Ojf ^2 



(1^1)) . 



(37) 



where A 



?( 



are the spatial eigenvalues of the operator A^. We use the following formal representation to eliminate 



the logarithm^ 



ln6= -lim / _e^*(^+^^). 



e^O /q t 



(39) 



Eq. (|37l) can be cast into the form 



I I d^k r+°° duj ^ '•+°° 



"t^ j.^^ ; dt -it(^^~\l{\k\)+ie^ _ i r d^k 



dt ^ 



— lim 



+ 00 



-»t(^-A^(|fc|)+»e) 



2tt 



iN r d^k 



Elim 



where we have used the following representations 



E 



d^k 



3 \/ i 



in) 



A 

8^' 



dw 
2^' 



N 



and 



-t-oo 



(40) 



(41) 



with Im (a;) > and z = — i. Note the presence of the redshift function in Ea. (|40p . This is a remnant of the 
original time component. In order to evaluate the integral over momenta in Eq. (j40p . we use the WKB approximation. 



^ This result is in agreement with the result of Ref. [Tsl , where only the graviton contribution contributes to the evaluation of the effective 
action. 

2 

ln-=lim/ ^^it{a+i.) _^it(b+u) (3g) 
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With the help of Eqs. (|B19l) . we define two r-dependent radial wave numbers ki (r, I, Ai^„/) and ^2 (r, I, X2.ni) for the 
Lichnerowicz operator (TT tensor) 



fc? (r,Z,Ai,„0 = A?^„,-^-m? (r) 

kUr,lA2,ni) = \U-^-ml (r) 
and we separate the effective masses in two pieces 

'^1 ('') = "^i (^) + "^1 5 (^) 



(42) 



with 



and 



m| (r) = m| (r) + mj g (r) 



6(r) 



mi (r) = 3 1 



mj s (r) 
"i2,5 (^) 



^6'(r)-^6(r) 



2r' 



' b'ir) + ^bir) 



(43) 



(44) 



(45) 



The WKB approximation we will use is equivalent to the scattering phase shift method and to the entropy computation 
in the brick wall model. We begin by counting the number of modes with frequency less than A^, i — 1,2. This is 
given approximately by 



(46) 



g(A,)= / ly,{l,X^){2l + l)dl, 
Jo 

where Vi (Z, A^) is the number of nodes in the mode with {I, Ai), such that (i — 1, 2) 

1 f+°° I 

Vi{l,\i)^-i dx\ kf {r,l,\i). {r = r{x)) 

J-00 ^ 



(47) 



In Eq. (j47p is understood that the integration with respect to x and / is taken over those values which satisfy 
k1 (r, Z, Xi) > 0, i = 1, 2. With the help of Eas. (|46l 47p . the total energy associated to the energy density in Eq. pO)) 
becomes (r = r [x)) 



— ^ / dxN (r) r2 / ^ A,2 ^ Xj - m| {r)dX,. 



By extracting the energy density, we obtain 



A 
8^ 



Xi\X'l-ml{r)d\, 



(48) 



(49) 



where we have included an additional 47r coming from the angular integration and where we have included in the 
volume term the redshift function. Of course, Eq. (j49p is divergent and must be regularized. 

IV. REGULARIZATION AND RENORMALIZATION OF ONE LOOP CONTRIBUTION TO THE 

COSMOLOGICAL CONSTANT 



We adopt the zeta function regularization scheme and by introducing the additional mass parameter ^ in order to 
restore the correct dimension for the regularized quantities, we define 



1 

4^' 



,2£ 



+00 



dXi 



A? 



1,2. 



(Af - mf {r)f 



(50) 
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The integration has to be meant in the range where Af — mf (r) > 0. FoUowing the same steps as in Ref.[i^, one gets 

In order to renormahze the divergent ZPE, we write 

J_ A(^ A*-^j^ mj {r)+mf (r) 

SttG 8nG SttG SttG diir^e ' ^ ' 

Thus, the renormahzation is performed via the absorption of the divergent part into the re-definition of the bare 
classical constant A. The remaining finite value for the cosmological constant reads 

To avoid the dependence on the arbitrary mass scale /i in Eq. (|53p , we adopt the renormahzation group equation and 
we impose that [1^1 

Solving it we find that the renormalized constant Aq should be treated as a running one in the sense that it varies 
provided that the scale fi is changing 

Aq (^,r) _ Aq (^o,r) ^ m\ {r) + m\ {r) ^ ^ 

SttG ~ SnG + 32^^ ^""jTo- ^^^> 

Substituting Eq.(I5S]) into Eq.(IS31) we find 

Potentially, we have three cases: 1) to| (r) 3> m| (r), 2) m| (r) — m| (r) and 3) to| (r) <^ m| (r). Case 2) reduces 
to a single point and therefore will be discarded in this analysis. In case 1) essentially we consider a long range 
contribution of the graviton which will be vanishing for r — > oo. Finally, case 3) is a short range case and it leads to 



2 

8ttG ^ 647r2 ^ 

Z— 1 



4 ( 

I, o Ir) in — 



(57) 



The above expression works for a background described by Eq.Q which must satisfy the Einstein's field equations. 
We specialize the result to the case of interest, namely the Schwarzschild and the de Sitter metrics. However, since 
the exterior part of the gravastar is of the Schwarzschild form, we begin with this background which is in common 
with a wormhole model. In the short range approximation we find 

"^?,s W = -^2,8 W = W = ^^3^ (58) 

and the range of validity is when r e [2A/G, 5MG/2] which can be determined by case 2). In this range, 
Aq {ho, r) I (SttG) has the following properties: 

i) For r = 2Af G, Ao (^io, 2MG) oo when M ^ 0. 

ii) For r > 2MG, Aq (pQ,r) when M -> 0. 

This can be summarized in the following double limit 

lim lim AninQ.r)^ lim lim AninQ.r) (59) 

which appears to be a sort of non-commutativity appearing in proximity of the throat. A similar behavior was 
conjectured by Ahluwalia[l3| in connection with the black hole entropy where a relation of the type 

[Ipls] = iXl (60) 

was introduced. The "Schwarzschild" Is and the "Planck" Ip lengths are no more simply lengths, but operators. An 
analogy can also be found for Yang-Mills theory in a constant chromomagnetic background [iSj, even if the situation 
interests the infrared region instead of the ultraviolet region. This non-commutative effect of Ea. (|59p could be 
interpreted as a signal of a phase transition. 
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A. Persisting of the wormhole 

With the help of relations ([58)) and from Eq. ([57)) . setting r = rt = 2AIG, we find that Ea. ((57| becomes 



' mf^l, (61) 



8ttG USiT^rf \8r^fil 

namely the throat does not manifest quantum fluctuations. Ea.(|5T|l satisfies the following inequality 

Ao (po.n) Aoino, ft) 9 /Se _ 

We recognize that the expression of the upper bound is of the Casimir form, in the sense that we have a computation 
procedure mimicking the Casimir device whose plates are located at the throat and at infinity. As in a Casimir device 
whose energy density is proportional to the inverse fourth power distance of the plates, also this case manifests the 
same behavior. From this it is evident that the gravitational field of the wormhole can never be switched off. However 
as M becomes smaller and smaller, one cannot avoid to enter in the quantum phase of the throat. Therefore, at a 
certain distance from the throat, a "brick wall" [l9| can be formed due to quantum fluctuations. This forbids the 
throat to be reached. In some sense, it is the Casimir energy that changes the structure of the throat. If this is the 
case, Eq. (|6T1) becomes 



(63) 



where rt is of the form rb = rt + h with h representing the "brick walF . The ZPE is now regular when the throat 
vanishes. 



B. Wormhole turning to a gravastar 

In this case, the Casimir energy not only creates a "brick wall", but changes completely the structure of the 
wormhole by a topology change^. The effect of the vacuum reorganization induces a redefinition of the function b (r) 
in © 



1 2MG(^^y 0<r<ro 
\ 2MG r>ro 



bir)^{ ^^"^V^y u^, ^,0 (64) 
The "inner cosmological constant" is of course 

A. - ^ (65) 

and regulates the small de Sitter universe inside the wormhole. Of course the new shape function is continuous in rg, 
but it leads to a completely new scenario, because it converts a wormhole throat in a "cosmological throat" , even if 
of very reduced size. The metric is now regular at the origin and even if b (r) is continuous at rg, the energy density 
is not. This is in agreement with the initial setting of Eq.([2]), but without the thin shell. The same discontinuity 
reappears to one loop level for in the inner region, Eq. (j57p becomes 



Ao(Aio,r-o)|„, 1 f6MG\\ f6MGy^\ 

SnG - -32^ [-^) [-^) ' ^''^ 



while in the outer region one gets 



SnG - -32^ [^) [-^) • 



^ See also Ref. [20l| where a discussion on the possible topology change induced by Casimir energy between dark stars and wormholes is 
faced. 
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However, if we try to insist and impose continuity on the boundary 

AoImcJ-o)!™ = Ao(Afo,»'o)|o„t, (68) 
we find a solution provided one fixes the renormahzation point to 

3Af GV^ 2 



= Mo- (69) 



C. Comparing the gravastar to the wormhole 

We are now ready to compare the gravastar with the wormhole model. As discussed in the introduction, we know 
that both models share the same ADM mass. Essentially, the ADM energy is defined as 

Eadm = ^ {D'h,, - D,h) r= , (70) 

where the indices i, j run over the three spatial dimensions and 

hi] = gij -9ij, (71) 

where gij is the background three-metric. Dj is the background covariant derivative and is the unit normal to the 
large sphere S. However, Hawking and Horowitz[2lj have shown that the definition ([70]) is equivalent to 

Eadm = ^ / d^x^ {k - fc") , (72) 

where a is the determinant of the unit 2-sphere. represents the trace of the extrinsic curvature corresponding to 
embedding in the two-dimensional boundary ^5' in three-dimensional Euclidean space at infinity. In Eq. (j72l) . it is well 
represented the subtraction procedure between two metrics having the same asymptotic behavior. Therefore with 
a natural extension we define the subtraction procedure in such a way that we can include quantum effects: this is 
the Casimir energy or in other terms, the vacuum energy. One can in general formally define the Casimir energy as 
follows 



[dM] = Eo [dM] - Ea [0] , (73) 

where Eq is the zero-point energy and dM. is a boundary. For zero temperature, the idea underlying the Casimir 
effect is to compare vacuum energies in two physical distinct configurations. The extension to quantum effects is 
straightforward 

Ecasimir [dM] = {E, [dM] - E, [0]),;^,,,,,; + {Eo [dM] - Eo [0])i_;„„^ + . . . . (74) 

In our picture, the classical part represented by the AVM-li^e energy is vanishing, because the asymptotic behavior 
is the same for both the wormhole and the gravastar. This means that 

Ecasimir [dM] - {Eo [dM] - Eo [0])i_;„<,p + ...., (75) 

namely Ecasimir is governed by purely quantum fiuctuations. Here with Ecasimir we mean that the energy is an 
energy density. Thus, the Casimir energy can be regarded as a measure of the topology change, in the sense that 
if Ecasimir IS positivc then the topology change will be suppressed, while if it is negative, it will be favored. It is 
important to remark that in most physical situations, the Casimir energy is negative. Consider now the one loop term 
and suppose to compare a gravastar and a wormhole with the same mass M and the same renormahzation point /ig. 
If we take flat space as a reference space, we can write 

« [dM] - E^^ [dM]) = « [dM] Eo [Q]),_,^^^ + {Eo [0] - E^^ [dM]) , (76) 

where Eo [O]i-ioop represents the ZPE contribution of fiat space to one loop, which is absorbed into the regularization 
procedure. Since outside the gravastar radius rp 

E^ [dM] = E^^ [dM] , (77) 
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we consider the region in proximity of tq. There remains to be evaluated the density energy difference 



-.^(A.(«--- 



loop 



,n) 



Ao {^J-o,rQ) 



(78) 



where Aq (/io, rb)^ means that we are evakiating the one loop term due to the wormhole background in proximity of 
the brick wall, where we expect to receive the largest energy density contribution. The minus sign appears because 
of a consequence of the definition of the induced cosmological constant of Ea. ((T2|) and the definition of the energy 
density of Eq.®. We have to discuss when 



E}^ \dM] - -Bn^* \dM 



^0 



GS 







' 1 — loop < 



^0. 



After algebraic manipulation, we find that this happens when 

.3 X 4 



In 



where we have used Ea. (j69p and we have defined 



■In 2 



>0 0<x<l; x > 1.26 
= x=l; 2; = 1.26, 
< 1 < X < 1.26 



(79) 



(80) 



n 



y 1 




FIG. 1: Plot of {E"^ \dM\ 



'dM\ 



-loop 



as a function of x = ro/rt- 



The situation is better illustrated in Fig.(IT]), where we immediately recognize that a tiny region exists where the 
ZPE is negative. This means that, in this range, the permanence of a wormhole is energetically favored with respect 
to a gravastar of the same mass M. The situation changes significantly if we avoid to fix the renormalization point 
with the choice of Eq. (|69p This means that we are abandoning the continuity between the external and the internal 
region of the gravastar to one loop. In this case, Ea. (l79l) becomes 



where x is given again by Eq. ()8ip and 



y 



41n2 - 3x^ In: 



3MGVe 



(82) 



(83) 



The behavior of the ZPE difference is shown in Fig. ([2]), where the equality of expression ([7^ is reached when 
y = X = 0.4948. It appears that below y and x, the ZPE becomes negative denoting that the ^^ground state" between 
a wormhole and a gravastar is represented by a wormhole whose radius has to be smaller than the gravastar radius. 
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-4- 



-5-" 

FIG. 2: Plot of {E^ [dM] - E^'^ {9M])^_^^^^ as a function of x = ro/rt and y. 



V. CONCLUSIONS 

In this paper we have considered one loop corrections to the AVAi mass of a wormhole and a gravastar respectively. 
The motivation comes from the fact that a gravastar can represent an alternative to a black hole or a wormhole as 
regards the gravitational collapse. Nevertheless various problems of comparison arise because outside the gravastar 
radius, the metric is of the Schwarzschild type. If we adopt the energy point of view we find that, asymptotically, 
these configurations share the same AVAi mass and therefore even in this case they are indistinguishable. However, 
since the gravastar has a different core with respect to the wormhole, a difference between them can emerge from a 
ZPE contribution. To this purpose we have computed an effective action to one loop, by looking at perturbations 
on the space-like hypersurface S. Since the perturbation involves only the spatial part of the metric, ghosts do not 
come into play for the energy contribution. Therefore, only the graviton is important to establish what happens to 
the ZPE. Note that every form of ZPE can be interpreted as an induced cosmological constant. In our case, this 
interpretation is very useful to apply standard regularization and renormalization procedures. In a sense, we can 
think that the ZPE induced by a wormhole or a gravastar contributes to a cosmological constant. On the other 
hand, we can think that the induced cosmological constant could be used to give a sorting to ZPE. It is interesting 
to note that, it is the double limit ((59|) that denotes that something different appears in the throat proximity. This 
difference is principally caused by ZPE or Casimir energy which becomes so intense to create a thick barrier ("brick 
wair) or a topology change (^'gravastar"). Therefore, it is important to discover under what condition a wormhole 
persists or change into a gravastar or vice versa. It appears that a fundamental element to understand in which 
direction the geometry becomes relevant is in the radii ratio (j8ip . Note that the comparison between the gravastar 
and the wormhole is done with the same mass M, as it should be. The other important parameter seems to be the 
renormalization point which, in this example, translates the presence or the absence of an energy gap between the 
external and the internal region of the gravastar. If one imposes the continuity of the ZPE through the gravastar 
radius, one meets a constraint on /iq leading to the plot of Fig.([T]) showing that the region of permanence of the 
wormhole is very subtle. Indeed, from inequality (|M|. outside the range 1 < a; < 1.26, the wormhole appears as an 
"excited state^^ with respect to the gravastar. In particular, one should note that the minimum for the ZPE difference 
appears for x — 1.152 with a minimum value of —0.244. On the other hand, if one abandons the continuity condition 
of the gravastar ZPE and treats /zq as a free parameter, one finds from Fig. ([2]) that the region of permanence of the 
wormhole is larger. Therefore we arrive at the conclusion that the permanence of a wormhole or a gravastar in their 
reciprocal comparison is strictly related to their size. 
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Appendix A: Disentangling the gauge modes 

To explicitly make calculations, we need an orthogonal decomposition for h^i, to disentangle gauge modes from 
physical deformations. To this purpose it is convenient to decompose h^i, into a trace, longitudinal and transverse- 
traceless part in D dimensions [H, [2^ [l^: 

h^. = ^hg,,, + (LO^, + (Al) 
where the operator L maps into symmetric tracefree tensors 

(^eU-V^e. + V.C^-|.g^.(V-e) (A2) 

and 

g'^'-h^, = 0, = 0. (A3) 

The decomposition (|A1[) is orthogonal with respect to the following inner product 

M 



{h,k):^ I d''x^G>'''P"h^Ax)kp,{x), (A4) 

where 



Q^^<^P<T ^ li^gf^Pg'^'^ + g^^'^g'^P + C g^" gP") (A5) 

and C is a constant. For the positivity of G^'^'"^ C > —1/2. The inverse metric is defined on cotangent space and it 
assumes the form 

(p, q) / d^'x^G.^p^p^^'' (x) (x) , (A6) 
Jm 

so that 

G^'^^^G^p,, = i {S^S:; + SiiS;) . (A7) 

Following Ref.jl^l, we observe that under the action of infinitesimal diffeomorphism generated by a vector field e^, 
the components of (|Aip transform as follows 

^M^^M + ^M' h^h + 2V^^e^, hj^u^h^.- (A8) 
We can fix the gauge freedom (lASp , by fixing 

Cm = 0. (A9) 

if the manifold admits conformal Killing vectors, namely vectors annihilated by the operator L, one additional gauge 
condition involving the trace is necessary. Assume that such vectors are absent form the manifold. The Jacobian 
factor induced by the change of variable, namely h^^ — >■ (/i, /i^^) is 

J = det(LtL)^, (AlO) 
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where the determinant is calculated in the space of vector fields excluding conformal Killing vectors. Using the 
definition (PO]) . the operator acting on vector fields is 



We can write decomposition (|Aip in the following way 



(All) 



(A12) 



The change of variables h ^ a does not introduce any additional Jacobian factor, then the path integral measure 
separates into 



Vh^,^ = det ipL) " Vh^ViVa 



and the quadratic part of the action (|T8|) becomes 

1 



2k 



M 



4 ^ 32 



(A13) 



(A14) 



The path integral can be written in terms of functional determinants of transverse-traceless tensor fields (indicated 
by T), vector fields (indicated by V) and scalar fields (indicated by S): 



Z = det (Ai)"3 det {L^L) ' det (-A)" 



(A15) 



The determinant over vector fields is the analogue of the Faddeev-Popov determinant. We can further decompose it 
by introducing a Hodge decomposition 



(A16) 



where -0 is a zero-form, is a harmonic one-form and a; is a two-form. Excluding the presence of harmonic vectors 
we find that the operator (jAlip separates into 



and using the equations of motion on shell R^^, = 0, we get for the Jacobian 

det (L^'LY =det(A)5det(A)^ 



VII 



and the one loop effective action simplifies 



^i-ioop = ^ [Trln Af-- Tr In Ayi] = ^ 



M 



d^k 



(27r 



— In Xj^j^ 



d^k 



In Xi 



(A17) 



(A18) 



(A19) 



where X^j. and Ay^ are the eigenvalues of the Lichnerowicz operator for TT tensors and the transverse vector operator 
respectively. Before going on, we need to precise a point concerning the computation of a functional determinant. 
Generally speaking the functional determinant of a given differential operator O can be represented by 



det O = exp [Tr In O) . 
However within the zeta function regularization, it is no longer true that 

det AS = det^ldetB, 
where A and B are two elliptic operators. In general, one has 

det AB = exp a {A, B) det AActB. 



(A20) 
(A21) 
(A22) 



where a (A, B) is a local functional called multiplicative anomaly. As pointed out in Ref. [2J| , one can assume the 
multiplicative anomaly to be trivial, namely a (A, B) = 0. This is justified by the fact that to one loop approximation 
a non-trivial multiplicative anomaly may be absorbed into the renormalization ambiguity. 
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1. Evaluating detv {L^L^ in 3 dimensions 



In Eas. (jAl[) . (jA2l) and (jA17p . we have simply to put D = 3. In addition, we have to remind that, in this case the 
Ricci tensor is not vanishing. If is further decomposed into a transverse part with V"^J — and a longitudinal 
part with ^| = Vqi/i, then the orthogonal decomposition reduces to 



(A23) 



M 



M 



(A24) 



with 



A'^^ = l\g^^ -R-K 



(A25) 



The decomposition (IA24[) is orthogonal up to the — ■(/'-mixed terms. If we try to compute the related Jacobian 
induced by the vector-scalar part in Ea. (|A24[ ). we obtain 



exp 



Ji j V^Vipcxp 
In terms of the functional determinant 



M 



= 1. 



Ji = 



(det^T^^, A/(*'^') 



(A26) 



(A27) 



where M^*'^' a (3 -|- 1) x (3 -I- l)-matrix differential operator whose first 3 columns act on transverse spin one field S^j 
and whose last columns acts on the spin zero field ?/'• From Eq. (jA24[ ). the corresponding matrix can be read off, 



-4i?'JV, 



-2 (I A2 +V,i?*J Vj) 



(A28) 



Thus Eq. (|A18p becomes 



(^detjT^^ 



M 



= (det A^) ' det 



^ A2 +V,ff^ V, 



(A29) 



The presence of the Ricci tensor in the scalar term of Eq. (jA29|) is an artifact of the foliation, because in 4 dimensions it 
disappears. By means of the contracted Bianchi identities V^i?*-' = R/2, we can simplify somewhat the expression 
of the scalar term of Eq. (jA29|) 



- A^ +ViR'^Wi = - A^ +ivJ'i?V, + R'^'^iV, 
3 3 2 

The determinant of the operator in Eq. (|A30p can be cast into the following form 



(A30) 



det ( ^ A^ +ff^ V,Vj j ' = exp ^Tr In A^ +i?*^ ViV. 



exp 2"^'" 



-A' 



-A' 



t V J 



exp -Tr 



In 1- 



-A^ 



(A31) 
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The second term is a correction to the principal part and will not be considered in the W.K.B. method. The main 
part of the operator reduces to 

det 0^2^ ' = (det A)5 , (A32) 

where we have absorbed the constant factor into the definition of the determinant and we have redefined the scalar 
wave function in such a way to absorb the operator A. To summarize, Eq. (|A18l) in 3 dimensions becomes 

det {L^Ly = (det A^)^ (det A)^ . (A33) 

With the help of the functional determinant in Eq. (|A33p , the one loop effective action can be written as 

Ti-ioop = (det A^) ^ det (A)^ (det A^,) (det A,) ~" . (A34) 

One can observe that the scalar determinants should cancel each other except for three subtleties. First the sign of 
the operator Act in (|30|) appears to be different from that in the Jacobian, second the integration leading to (IA29[) 
excludes zero modes not included in the Jacobian, so any cancellation will not be complete. This also happens in the 
full covariant computation. Last but not least, the cancellation should be done after integration over the time part 
in the determinant of the operator. 

Appendix B: The Lichnerowicz operator for TT tensors 

Our starting point is the expression (fT6|) and the metric For the benefit of the reader, we recall the representation 
of the operator O 

and we simplify the expression of the Riemann tensor in 3 dimensions 

R 

Rikji = gijRki - guRkj - gkjRii + guRi] - {gi^gu - gugkj) ■ (B2) 

Then, the operator O**^.?' becomes 

-VVa.g'^'g^' + -^9'^9^^ ~ 2 {g^m^^ - g^R''^ - g^m^^ + g^^W^) + R {g^^ g^^ - g^^g^^) 

^K^g^^ + W^g"^ - AR'^g''^ + Rg'''g^' = -V^VaO^'^^o^'' + ^g''"g^^ 

Ot'^ 

~ 2 {g'^R''^ + g'^'R'^) + R [g'^ g^^ ~ g'^g^^) + g^^ + 3i?^'^/ - g^^ + Rg^'^g^K (B3) 
When we fix our attention on TT tensors, we obtain a further reduction 

O'""'' = -V^VaS'^ff^' + + 2R''g^'. (B4) 



Thus the related eigenvalue equation is 



^^^^)\ + ^^H->?K (B5) 



where 

{^2h'^)\ := (-At/i^^)J + 2 (i?/i^^)J (B6) 
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and 



- (At/.-): = -As {h-i: + ^ (i - v^) i^n: ■ 

As is the scalar curved Laplacian computed on the background of metric ([3]), whose form is 



A. = 1 - ^ ^ 



6(r)\ (f fAr-b'{r)r-3b{r)\ d L 



2r2 



r I dr^ 

and i?° is the mixed Ricci tensor whose components are: 

'b'{r) b{r) b'{r) b{r) b' (r) b{r) 



dr 



Rf 



(B7) 



(B8) 



(B9) 



We will follow Regge and Wheeler in analyzing the equation as modes of definite frequency, angular momentum 
and parity psjl- In particular, our choice for the three-dimensional gravitational perturbation is represented by its 
even-parity form 



/i?^™ {t, r, d, <P) = diag 



H{r) 1- 



b{r) 



, r^K (r) , sin^ -dK (r) 



Yi^ {d, cj)) F (t) . 



(BIO) 



For a generic value of the angular momentum L representation (|B10[) . together with Eas. (jB7l BSP leads to the following 
system of PDE's 



-A, 



; + 2 (^^IW - ^) + ^) i/ (r) F it) = Xl,H (r) F (t) 
(- A, +2{^+'4r^)+^)Kir)F it) = Xl,K (r) F it) 



where A; is 



A, = 1 - 



6(r)\ d 



dr'^ 



The action of A; on the reduced fields 



/4r-6'(r)r-36(r)\ d lil + l) 6 / b{r) 
\ / dr r 



(Bll) 



(B12) 



IS 



A, 



/i,2 (r) 



Hir) 



Mr) 



Kir) = 



h{r) 



bir)\ £_ 
dr^ 



-2 1 



bir)\ f Ar - b' ir) r - 3b ir) 
~r) ^ V "ir 



d 
dr 



(B13) 



+- 



2 1 



bir) 



4r - b' (r) r - 36 (r) 
Yr 



/(/ + 1) 6 



bir) 



ha ir) 



1 



1 



dr"^ r 



b' ir) b ir) 



d 1 

dr r^ 



2 2r 

and using the proper geodesic distance from the throat 

dx = ± 



b' ir) b ir) 
~ 2^ 



dr 



U^ + 1) 6 / bir) 



b{r) 



hair) (B14) 



(B15) 
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we get 



1 



h' (r) h [r) 



2 2r 

where r = r {x). Thus, the system (|B11[) becomes 



I {I 



il _ 1 f 1 _ Mi;) 



ha (r) 



with 



^ + V,{r)\h{x)F{t)=\l,h{x)F{t) 
'£, + §^ + V2 (r)l h (x) F (t) = Xl j2 ix) F (t) 



Vi (r) = ^ + C/i (r) 
V2 (r) = ^ + C/2 (r) 



and 









-^b{r)] 


= m\ (r) 


[ U2 (r) = ^ 






^i^b{r)] 


= ml (r) 



where we have defined two effective masses dependent on r. Now let us consider the following state 



1 



{t, r, 9, <p\uj, k.l.m) = -fi {x) Yi^ {■&, <p) exp {-iujt) , 



then the system (jB17p becomes 



+ ^1 {r)\ h (x) = XlJi (x) 
+ V2 (r)l h (x) = Xlj2 (x) 



(B16) 



(B17) 



(B18) 



(B19) 



(B20) 



(B21) 
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